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Flowfield measurements are presented for a symmetrical NACA 64A010 airfoil section at transonic con- 
ditions. Measurements were obtained for three angles of attack with the freestream Mach number fixed at 0.8. 
The cases studied included a weak shock-wave/boundary-layer interaction, an interaction of medium strength 
with mild separation, and an interaction of sufficient strength to produce a shock-induced stall situation. Two 
nonintrusive optical techniques, laser velocimetry and holographic interferometry, were used to characterize the 
flows. The results include Mach number contours and flow angle distributions in the inviscid flow regions, and 
turbulent flow properties, including the turbulent Reynolds stresses, of the upper surface viscous layers, and of 
the near-wake. The turbulent flow measurements reveal that the turbulence fluctuations attain equilibrium with 
the local mean flow much faster than previously expected. 


Nomenclature 

c = chordlength of airfoil 

C L = lift coefficient, L{ 

C p = local pressure coefficient, (p— ) IVip^u 2 ^ 

Kq _ Q = Gladestone-Dale constant 
f = Prandtl mixing length 

L eff = effective optical path length 

M - Mach number 

N = total number of velocity realizations 

p = pressure 

q = turbulence kinetic energy, ((u’ ) 2 + (v' ) 2 + 

<h>') 2 )/ 2 

Re c - Reynolds number based o n airf oil chordlength 
R uv = uv correlation coefficient, u'v'/(u'}( v') 

u - velocity component in streamwise direction 
v = velocity component in normal direction 

w = velocity component in cross-stream direction 
x = coordinate in the streamwise direction 

y = coordinate in the normal direction 

a - angle of attack 

6 - viscous layer thickness 

e = fringe shift in interferogram 

0 = flow angle, arctan v/u 

X = laser wavelength 

p - fluid density 

r = Reynolds shear stress 

co = weighting factor for velocity biasing 

Subscripts 

t,2 = two arbitrary points in the flow 

oo = freestream conditions 

e = conditions at edge of viscous layer 

/ = rth velocity realization 

s - surface 

t = stagnation condition 

w = conditions at the surface 
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Superscripts 

( )' = fluctuating quantity 

( ) = time-averaged quantity 

<') = rms value of quantity 

Introduction 

A LTHOUGH considerable advances have been made 
toward the numerical prediction of the transonic flow 
past two-dimensional airfoil sections (e.g., Refs. 1-5), these 
numerical methods currently have a limited range of ap- 
plicability, either because they ignore viscous effects 
altogether or because they inadequately predict the viscous 
effect as the shock wave strengthens on the airfoil’s upper 
surface. Accurate predictions are especially difficult to obtain 
once shock-induced separation occurs. This inability to ac- 
count accurately for viscous effects is a direct consequence of 
the deficiencies in the models employed for the turbulent 
Reynolds stresses (e.g., Ref. 4). To assess truly the ability of 
these methods to describe the flow behavior and to provide 
insight into how to improve on these methods, quantitative 
measurements of the external flow and, especially, of the 
turbulent transport properties, are needed. In transonic 
testing, flowfield data are also needed if account is to be taken 
of wall interference effects which can alter significantly the 
pressure distributions from those that would be measured 
under free-air conditions. 

The objective of the present study was to obtain detailed 
descriptions of both the inviscid and viscous flow regions 
including the turbulence transport properties about a 
geometrically simple, symmetrical airfoil section (NACA 
64A010) for several different test conditions. The conditions 
selected for study were: a weak shock-wave/boundary-layer 
interaction for which inviscid methods should be adequate, an 
interaction of medium strength containing some separation, 
and a strong interaction with massive separation as occurs at 
stall. The severity of the interaction was established by 
varying angle of attack at a fixed freestream Mach number of 
0.8. To describe the flow about the airfoil section quan- 
titatively, two nonintrusive optical methods were employed. 
A laser velocimeter technique was used to provide detailed 
localized information of the streamwise and normal velocity 
components (mean and fluctuating). Holographic in- 
terferometry was used to obtain the density fields about the 
airfoil section from which Mach number contours were 
determined. For completeness of the experiment, surface 
pressure distributions, wake pressure, oil flow visualizations, 
and conventional Schlieren and shadowgraphs were also 
obtained. 
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Experimental Technique 

The experiments were conducted in the Ames 2- by 2-Foot 
Transonic Wind Tunnel which is a closed-return, variable- 
density tunnel with 21 % open porous-slotted upper and lower 
walls for transonic testing. A 6-in. chord 64A010 airfoil which 
spanned the entire test section was selected for the tests. 
Tunnel conditions were fixed at a freestream Mach number of 
0.8 and a chord Reynolds number of 2 x 10 6 . To insure that 
the boundary layer was turbulent at the foot of the shock 
wave for all test conditions, transition strips were affixed to 
the airfoil section at the 17% chord station on the upper and 
lower surfaces. The 1 .25-mm-w'ide strips consisted of 0.13- 
mm nominal diameter glass bends, selected in accordance 
with the recommendations of Ref. 6. Transition effectiveness 
was verified by the sublimination technique. 


Laser Velocimeter 

Thejaser velocimeter was a dual-color system utilizing the 
4880 A and 5145 A lines of an argon-ion laser; one spectral 
line to measure the streamwise velocity component, the other 
to measure the vertical velocity component. Bragg-cell 
frequency shifting, necessary for probing highly turbulent and 
separated flow regions, was incorporated in both spectral 
lines. The frequency offsets also facilitated the direct 
measurement of the vertical velocity component (i.e., ±45 
deg beam orientations to resolve the vertical velocity were 
unnecessary). The effective sensing volume was ap- 
proximately a cylinder, 200 pm in diameter and 3 mm long, 
w'hose axis was aligned with the cross-stream direction. 

Signal processing was accomplished with single-particle 
burst counters and the individual realizations from the two 
channels were recorded simultaneously with a digital com- 
puter. This allowed the velocity correlation mV to be ob- 
tained in a straightforward manner by multiplying and 
averaging, rather than by the less accurate method of sub- 
tracting the variances of the two signals obtained from ±45 
deg beam orientations. In addition, it facilitated the weighting 
of the data for the velocity biasing which is believed to occur 
when the particle concentration is low-. 7 This biasing is at- 
tributed to particle arrival rate dependency on instantaneous 
velocity. In this paper, the data presented were adjusted using 
the two-dimensional weighting factor suggested in Ref. 7. The 
mean velo cities , turbulence intensities, and velocity 
correlations wV w'ere calculated as follows: 


<V .V 

£ OljuJ ^ CO, 
i = l 1 /=/ 


( t »,»;/ } 

c».) 


i = l 1 r- 

= / 7 



- uv 


uV = ( £ 0),U,V, £ <o,)-* 

with the two-dimensional weighting factor given by 


( 1 ) 

( 2 ) 

(3) 


co, = J/yfuj+ vj (4) 

The equations for 0 and < v ' > are identical to Eqs. (1) and (2), 
except Vi is used in place of u r When no correction for 
velocity biasing is applied, co, is taken to be unity; the data 
were also reduced in this manner. A comparison of the two 
sets of results showed that the same overall conclusions 
regarding the flow behavior would be drawn using either 
approach. 

Naturally occurring particles in the tunnel were used for 
light scattering. In this facility, lubrication oil within the drive 
system vaporizes and then later condenses in the tunnel circuit 
to provide a generous supply of scattering centers. Previous 
measurements 8 across a normal shock have shown that these 


particles are small enough in size (estimated to be 1 pm) to 
give very good response to a step change in velocity at sonic 
speeds. At each point in the flow, at least several thousand 
velocity realizations were used to calculate the flow properties 
given in Eqs. (1-3). Measured velocities at the edge of the 
viscous layers were generally repeatable to better than 1 % and 
agreed well with edge velocities inferred from the measured 
surface pressures. The primary cause for these variations was 
changes in the tunnel operating temperature. Under constant 
test conditions, the inviscid velocities repeated to within 
0.25%. In the viscous layers, mean velocities were repeatable 
to better than 2% of the edge velocity— including the ex- 
tremely turbulent separated regions. 

Holographic Interferometer 

The holography system was designed to utilize the existing 
tunnel Schlieren mirrors. The setup is based on the 
arrangement suggested in Ref. 9 for converting large wind 
tunnel Schlieren systems to holographic visualization systems. 
In the present investigation, the dual-plate method 10 was used 
wherein an exposure is made on a photographic plate under 
no-flow conditions and subsequent plates are made at the test 
conditions. After processing, the no-flow plate and one of the 
plates taken at test conditions are positioned in a recon- 
struction plate holder, illuminated with the reference beam 
and aligned for infinite-fringe interferograms. Vibrations in 
the system are relatively insignificant since the pulse duration 
of the ruby laser is extremely short ( = 20 ns) and motion 
between exposures can be eradicated in the reconstruction 
process. 

To reduce the infinite-fringe interferograms to density 
contours, tw-o-dimensional flow' was assumed in the 
relationship 


P2~Pl : 


^cff Kq - D 


(5) 


where p ; and p 2 are the densities at two locations, t is the 
fringe shift between the locations, L eff is the effective optical 
path length taking into account the sidewall boundary layers, 
X the laser wavelength, and A' g _d the Gladestone-Dale 
constant. Since the correction in optical path length due to the 
sidewall boundary layers was only 2%, the geometrical width 
of the tunnel could have been used without any significant 
loss in accuracy. With the interferometer aligned in the in- 
finite-fringe mode, each fringe corresponds to a constant 
density contour. Once the density at one point in the How is 
established, the remaining contours can be determined from 
Eq. (5). 

Due to the relatively large span of the tunnel, the sensitivity 
of the system to density changes is quite good for two- 
dimensional testing. At for example, one fringe 

shift corresponds to only about a 0.5% change in density and 
a corresponding change in Mach number of about 1%. This 
and the high aspect ratio (4) of the airfoil tested, which 
produced good two-dimensionality in the flow, contributed to 
the quality of data obtainable with this optical method. 

Results and Discussion 

Shown in Fig. 1 are the lifting characteristics (C L vs a) of 
the 64A010 airfoil section for M oe =0.8 as obtained from 
integration of measured surface pressure distributions. Two 
sets of data are shown in Fig. 1 : those obtained in the present 
investigation and those from a previous study 11 in the same 
facility. The slightly higher lift coefficients observed in the 
present study are believed to be a result of transition strip 
placement. In the study of Ref. 11, the transition strips w'ere 
placed nearer to the leading edge (6.1% chordwise station). 
This would produce a thicker turbulent boundary layer at the 
foot of the shock wave and, hence, a reduction in lift. Other 
64A010 data, such as those of Stivers, 12 are not included, 
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Fig. 1 Lifting characteristics of 64A010 airfoil section for = 0.8. 


NACA 64A010 
= 0 8 
o = 6.2" 

O PRESENT STUDY, R« c = 2(10) 6 



Fig. 2 Reynolds number effects in present study. 


since the intent of Fig. 1 is only to show the performance of 
this airfoil in this particular facility with boundary-layer 
transition imposed. At the subject Reynolds numbers, sub- 
stantially different results would have been obtained if 
natural transition had been employed, as in Ref. 12. It is 
Important to note from Fig. 1 the sudden decrease in the lift 
curve slope at an angle of attack of about 4 deg which strongly 
suggests the occurrence of a shock-stall situation. The three 
angles of attack concentrated on in the present study were set 
angles of 0, 3.5, and 6.2 deg. As is apparent from Fig. 1, the 
two lifting cases studied occur just prior to and after the 
sudden decrease in the lift curve slope. 

As seen in Fig. 1, the variation in C L due to Reynolds 
number was observed to be quite small in the study of Ref. 1 1 . 
To verify this, surface pressures were obtained for a = 6.2 deg 
and Re c - 4x 10 6 in the present study. These data are com- 
pared to data obtained at Re c - 2x 10 6 in Fig. 2. Included in 
Fig. 2 are results from Ref. 1 1 for Re c = 4x 10 6 , which show 
the effect of transition strip placement on the surface pressure 
distribution. It is evident that the surface pressure 
distributions for the two Reynolds numbers differ only 
slightly. The insensitivity to Reynolds number in this study 
should not be construed to imply that the data would compare 
well with full-scale Reynolds number results. In the present 
study, and in that of Ref. 11, the transition strip placement 
had the effect of nullifying the effects of natural boundary- 



Fig. 3 Infinite-fringe inlerferograms. a) a = 0 deg; b) cr = 3.5 deg; c) 
a = 6.2 deg. 


layer transition, thereby reducing the Reynolds number 
sensitivity. Obviously, the boundary-layer thicknesses relative 
to the chord of the airfoil in this experiment were much 
greater than for a flight situation where natural transition 
occurs. However, the objective of the study was not to obtain 
data valid for flight conditions but rather to obtain data that 
could impove our understanding of the basic flow phenomena 
and, consequently, our predictive capabilities. From this 
point of view, it is highly desirable not to have the flow be 
sensitive to the complex transition phenomena. 

Infinite-fringe interferograms obtained with the 
holographic interferometer for the three angles of attack of 
interest are presented in Fig. 3. As was discussed in the 
previous section, each fringe corresponds to a line of constant 
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Fig. 4 Surface pressure comparison between directly measured 
pressures and pressures inferred from interferograms, a = 6.2 deg. 


density. In the interferograms, the shock wave on the wing’s 
upper surface (upper and lower for the 0 deg case), the tur- 
bulent wake and the turbulent boundary layer downstream of 
the shock, as well as the stagnation region near the leading 
edge are readily apparent. The major contrast for the three 
angles of attack is the thickness of the boundary layer on the 
airfoil’s upper surface. For the two lifting cases (a = 3.5 and 
6.2 deg), the boundary-layer thickness on the upper surface at 
the trailing edge is, respectively, 2 and 6 times that for the 0 
deg case. Also apparent is the very dominant forward cur- 
vature in the shock wave near the airfoil for a = 6.2 deg, 
caused by the rapid thickening of the turbulent boundary 
layer. In the 3.5 deg case, oil flow visualization confirmed 
that the boundary layer separates at the foot of the shock 
wave (this was not the case for a = 0 deg); however, the shock 
curvature is very much subdued as compared to the 6.2 deg 
case. Bifurcation (two shock waves) was not observed in any 
of the cases studied. 

Although the interferograms are of value from a flow- 
visualization point of view, their greatest contribution comes 
from the quantitative density information they can provide. 
Ip principle, if the fluid density can be established at one point 
in the flow, the density can then be determined at any other 
point in the flow from Eq (5), provided the flow is two- 
dimensional. The validity of Eq. (5) in the interpretation of 
these interferograms was established by comparing surface 
pressure inferred from the interferograms with the directly 
measured centerline surface pressures. The inferred surface 
pressures were obtained by assuming isentropic flow in the 
inviscid region and constant pressure across the viscous layer. 
The laser velocimeter was used to establish the required 
density at one point in the flow. In Fig. 4, this comparison is 
shown for a = 6.2 deg. In the 6.2 deg case, the total pressure 
loss across the shock wave was not neglibile. A correction for 
this loss in total pressure as determined from the maximum 
Mach number immediately upstream of the shock wave and 
the shock-wave inclination is included in Fig. 4. This 
correction produced better agreement at the trailing edge. 

Inviscid Flow Properties 

Under the conditions of isentropic flow, lines of constant 
density are also lines of constant Mach number. In Fig. 5, 



Mach number contour plots obtained from the interferograms 
under this assumption are presented for the two lifting cases. 
Although there are total pressure losses across the shock 
waves, the effect in computing a Mach number is slight. For 
instance, for the most severe case, a = 6.2 deg, the maximum 
error in Mach number caused by neglecting total pressure 
losses is only AM = 0.01 . 

As seen from Fig. 5, only at the highest angle of attack was 
a very definite local region of supersonic flow observed 
downstream of the shock wave. Downstream of this sonic 
region the flow slowly decelerates. For all three conditions, 
the shock waves were of the weak solution. 

An inviscid flow property not available from in- 
terferometry is local flow angle 6 = tan ~ 1 ( v/u ) . For the three 
test conditions, horizontal scans were made with the laser 
velocimeter at Ayfc- 0.083 intervals above the airfoil. The 
vertical distance was referenced from the chordline at the 25% 
chordwise station, which was the point of rotation of the 
model. Shown in Fig. 6 are near- and far-fietd results for a = 0 
deg. Included are the flow angles predicted by the method of 
Ref. 3 using free-air boundary conditions. For the near field 
(y/c = 0. 167), the data show a larger flow angle change than 
that predicted by theory. The reason for the slight differences 
in the near field may be due to the small disturbance ap- 
proximations used in the method of Ref. 3. If the particles 
were unable to track the flow, too small a change in flow 
angle would be the expected result. In the far field, the ex- 
perimental data show larger angles in an absolute sense above 
the trailing edge than above the forward portion of the airfoil. 
This characteristic was even more pronounced for the lifting 
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Fig. 6 Inviscid flow angles. 


cases, as seen in Fig. 6. Inviscid theory does not predict this 
trend even when an attempt is made to model the tunnel walls 
a c is possible with the method of Ref. 3. Surprisingly, the 
measured far-field flow angles aft of midchord were not 
measurably different for the two lifting cases. To what degree 
this far-field flow angle trend is a result of viscous effects 
(possibly due to entrainment) or a result of wall effects needs 
to be investigated. Scans farther downstream, which would 
have been of interest, were not possible due to the limitations 
of the tunnel window size. Included in Fig. 6 are the near-field 
flow angles for a = 6.2 deg. Noteworthy is the flow direction 
change at the shock wave, which is consistent with the shock 
orientation and strength at that location. These results 
demonstrate the degree to which these two optical techniques 
used in conjunction can quantify the inviscid flow. 

Turbulent Flow Properties 

As a first step in exploring the viscous flow region with the 
laser velocimeter, a streamwise station just downstream of the 
trailing edge (x/c= 1.02) was surveyed. This station was 
sufficiently close to the trailing edge for the measured tur- 
bulent flow properties to represent the conditions at the 
trailing edge, yet sufficiently downstream that possible 
measurement difficulties associated with flare from the 
model's surface could be avoided. Figure 7 shows the mean 
velocity profiles obtained at this streamw'ise station for the 
three angles of attack selected for the study. The vertical 
distance is referenced to the trailing-edge location. Evident 
from this figure is the dramatic increase in boundary-layer 
thickness that occurs when the angle of attack is increased 
from 3.5 to 6.2 deg. (Recall that shock-induced separation has 
occurred in both cases.) Notice the upward displacement of 
the wake for the 6.2 deg case. This rapid displacement of the 
wake was also apparent in the interferogram. Surprisingly, at 
this streamwise station wake closure has already occurred, 
even for the a = 6.2 deg case. 

The corresponding turb ulenc e intensities <w' > and < v ') 
and the velocity correlation u'v' plotted against y/b are given 
in Figs. 8 and 9. The reader should note that each of the 



Fig. 7 Trailing-edge mean velocity profiles. 


velocity correlation profiles (Fig. 9) is plotted on a different 
scale. From Fig. 8, it is seen that the fluctuation levels in both 
the streamwise and vertical direction increase with increasing 
angle of attack; however, the increases are greater for the 
streamwise fluctuations. The largest changes occur in the 
central portion of the layers where the increase in anistropy of 
the flow with angle of attack is apparent. For all three con- 
ditions, the maximum fluctuation levels occur away from the 
wall, as would be expected for a shear flow that has ex- 
perienced an adverse pressure gradient. 

The Reyn olds shear stress for compressible f lows is given 
by — pu * v ' (assuming that the triple correlation, p ' u ' v ' , is 
negligible). However, in the transonic flow under in- 
vestigation, the change in p across the viscous layer is quite 
smal l (e.g ., p e fp w = 0.9 for M e = 0. 8) so the distribution 
of -u' v' /uj is nearly the same as - pu' v' / p w ui . Thus, 

-u'v' fu 2 ^ can, for practical purposes, be interpreted as the 
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M* = 0.8 
Re c - 2(10) 
x/c - 1.02 


6 


valid in describing the shear behavior. In Fig. 10, Prandtl’s 
mixing length f scaled to the boundary layer thickness 


a. deg 
O 0 



Fig. 10 Trailing-edge mixing lengths. 

nondimensional Reynolds shear stress. In the region where the 
flow is most energetic. Fig. 9 shows that the velocity 
correlation and, accordingly, the Reynolds shear stress, in- 
crease dramatically (by over an order of magnitude between 
the 0 and 6.2 deg cases), as do the turbulence intensities. This 
increase is very nearly equal to the change in turbulence 
kinetic energy, q 2 = ((u' ) 2 + (v' ) 2 + (w' ) 2 ) /2, where 
< w' > has been taken equal to <u' ). 

From Figs. 7 and 9 it is seen that the point of maximum 
shear stress occurs where the strain rate du/dy is the greatest 
This would imply that an eddy viscosity description may be 


( _/u'v' y idiu/u^) 

6 \ ui ) I d(y/&) (6) 

in the middle portion of the layers is compared to the model 
suggested by Escudier. 13 Since accurate gradients are ex- 
tremely difficult to obtain from experimental data, especially 
when the data are sparse, the determination of t was only 
attempted in the region where the slope was nearly constant 
(i.e., where d 2 ii/dy 2 was small) and thus, best defined In 
fact, one slope was established for each of the three cases and 
the mixing lengths were calculated only for those points for 
which this constant slope appeared valid. The behavior of f/6 
for ^/6<0.2 is not discernible from these data due to the 
uncertainty of the strain rate in this region. It is astonishing 
that the calculated mixing lengths agree so well with the value 
of 0.09, suggested by Escudier, in the middle portion of these 
layers. These three boundary layers have had entirely dif- 
ferent flow histories, yet in the central portion of these layers 
equilibrium between the mean flow and the turbulent fluc- 
tuations appears to prevail at the trailing edge. The boundary- 
layer thickness just upstream of the shock wave is extremely 
thin. Thus, the trailing edge is many initial-boundary-layer 
thicknesses downstream of the interaction region, which may 
account for the unexpected equilibrium between the mean 
flow and the turbulent shear. 

Measurements at several different streamwise locations 
were obtained for a = 6.2 deg. This angle of attack was 
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Fig. 12 Turbulence intensity distri- 
butions: a = 6.2 deg, =0.8, and 
Re c = 2 x 10 6 . 


selected for studying the streamwise variation of the flow on 
the basis of: 1) the thick viscous layer present, which 
facilitated detailed probing, and 2) the inability of current 
numerical methods to predict even closely a flow with such a 
strong inviscid-viscous interaction. These data are presented 
in Figs. 1 M3. For the measurement stations above the airfoil, 
the vertical distance is that from the surface, whereas in the 
wake it is referenced to the trailing-edge location. Along the 
airfoil, measurements were obtained within 0.5 mm of the 
surface. In Fig. 1 1 , the profile plots are vertically displaced in 
accordance with the airfoil surface location. The furthermost 
upstream measurement station was 0.30 c downstream of the 
separation point determined from oil flow visualization. 
Noteworthy from the mean velocity profiles are 1) the 
relatively small reversed velocities observed in the separated 
flow' region, 2) the free-shear-layer shape to the profiles at the 


stations along the wing’s surface, 3) the rapid wake closure at 
the trailing edge, and 4) the sudden upward displacement of 
the near wake. The small reversed velocities imply that the 
skin friction is near zero in that part of the separated region. 
The experimental evidence suggests that the flow is not wall- 
dominated, once separation has occurred. Beyond separation, 
the viscous layer appears to take on a free-shear-layer form 
with the outer edge of the layer following nearly a horizontal 
path. 

In Fig. 12, the corresponding streamwise and vertical 
fluctuation intensities are presented. At the most upstream 
station, the streamwise fluctuations, <w , )/w (X> , are seen to 
exceed the vertical fluctuations by more than a factor of 2. 
The streamwise fluctuations attain a level of over 40% at 
y/b = 0.5, indicating the presence of instantaneous negative 
velocities at an appreciable distance from the surface. 
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Fig. 14 Streamwise decay in production rale, turbulence kinetic 

energy, and shear. 

• 

Proceeding downstream, the streamwise fluctuations and, 
correspondingly, the turbulence kinetic energies are decaying, 
whereas the vertical fluctuations remain nearly constant. In 
the near wake (i.e., the last two downstream measurement 
stations), the character of the turbulence intensity 
distributions is seen to change rapidly; at the last station the 
flow has attained a nearly isotropic condition^ 

As seen from Fig. 13, the peak shear (-u'v'/ul is taken to 
represent pu'v'/p w ul) decreases with downstream distance 
as does the turbulence kinetic energy. At the last three 
stations, the presence of the lower surface boundary layer 
becomes evident from the sign change in the shear. Notice that 
the shear does not quite pass through zero where du/dy = 0, 
(i.e., w min location) indicating nonequilibrium effects. The 
lo catio n of the peak shear r maK (i.e., the location where 
-u'v'/ul, is a maximum) was found in all cases to 
correspond to the location where <w')/w» was also a 
maximum (see Fig. 12). 



Fig. 15 Distribution of uv correlation coefficient, R uv . 


A better illustration of the decay in the maximum shear 
stress and turbulence kinetic energy (O') taken equal to 
<u'» with streamwise distance is pre sented in Fig. 14. In- 
cluded in this figure is the quantity w' v' (du/dy) scaled to 
and c. This term is the primary source term of turbulence 
kinetic energy for this flow. As seen from Fig. 14, the 
production rate of turbulence kinetic energy decreases rapidly 
with streamwise distance. This is due in a large part to the 
decrease in strain rate, dti/dy , with increasing streamwise 
distance. With this rapid decrease in production rate, the 
corresponding rapid decay in turbulence kinetic energy is not 
surprising. These results, however, are contrary to those 
reported for a biconvex airfoil section at 0 deg angle of attack, 
where the shear and turbulence kinetic energy downstream of 
the shock wave and in the near wake were observed to increase 
with streamwise distance although the strain rate was rapidly 
decreasing. 14 The observations of Ref. 14 may be a result of 
flow unsteadiness caused by having separation on both lower 
and upper surfaces. 

The uv correlation coefficients R uv for the three upstream 
stations are given in Fig. 15. In this figure the zero-pressure 
gradient boundary layer results of Klebanoff 15 and the free- 
shear layer results of Leipmann and Laufer 16 are also shown 
for comparison. As seen from Fig. 15, the present data could 
be represented quite well by a single curve. In the inner half of 
the boundary layer the correlation decreases. This behavior 
has also been observed for incompressible boundary layers 
subjected to an adverse pressure gradient. 17 
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The mixing lengths at the points of maximum shear (scaled 
to 6) (this also corresponds to where the strain rate is a 
maximum) are presented in Fig. 16 as a function of distance 
along the airfoil. As shown earlier for the trailing-edge 
station, C/6 for the two upstream stations is well described by 
the value of 0.09 suggested by Escudier. In the wake, the 
mixing length scaled to 6 (the distance from the w min location 
to the upper edge of the wake in this case) is increasing. This 
behavior in C/6 is reasonable considering that for an in- 
compressible far-wake (76 = 0.18 (Ref. 18). The mixing 
lengths in the lower part of the wake where ItI attains 
another relative maximum are smaller as a result of the 
thinner lower surface boundary layer. These mixing lengths at 
x/c = 1.17 and 1 .33 are approximately one-third and one-half, 
respectively, those in the upper part of the wake, which in- 
dicates an even more rapid increase in mixing length with 
streamwise distance in this part of the wake. 

The Escudier model has been shown to be reasonably ac- 
curate in relating the Reynolds shear stress to the strain rate 
du/dy over the aft one-third of the airfoil’s upper surface for 
an extremely strong interaction case and at the trailing edge 
for weaker interactions. Although measurements were ob- 
tained only at the trailing edge for the two weaker in- 
teractions, the results for or = 6.2 deg strongly suggest the 
same would be true for these two cases. Yet, predictions for a 
64 A0 10 airfoil section using this turbulence model have been 
found to be inaccurate for lifting cases. 5 It appears that this 
mixing length model fails in the immediate vicinity of the 
shock wave. Extrapolation of the turbulence kinetic energy 
forward of the upstream-most measurement station in Fig. 14 
suggests extremely high fluctuation levels at the stagnation 
point. An accurate accounting of the Reynolds normal 
stresses at and near the separation point may be what is 
needed to predict the initial perturbations to the flow. Once 
the initial perturbation to the boundary layer is accurately 
modeled, the present results suggest that a zero equation 
model such as Escudier’s would be appropriate in the 
downstream regions if modifications were made to account 
for the increasing length scales in the nearwake region. 

Concluding Remarks 

Measurements of the flow behavior about a 64A010 airfoil 
section at transonic conditions have been presented. The 
results include Mach number contour plots and flow angles 
for the inviscid flow regions, and the turbulent flow 
properties of the upper surface viscous layer on the aft 
portion of the airfoil and of the near wake. 

In view of the apparent inadequacies of current models for 
the turbulent Reynolds stresses, the turbulent flow 
measurements presented must be considered the most relevant 
results of the study. The present results show that equilibrium 
between the turbulent fluctuations and the local mean flow 
occurs much more rapidly in the central portions of the 
viscous layers than expected with //5 = 0.9 in this region. This 
condition was found to hold even at the furthermost upstream 


measurement station (x/c = 0.61) for a = 6.2 deg. Although 
measurements were only obtained at the trailing edge for the 
two weaker interactions, the results for a = 6.2 deg strongly 
suggest the same would be true for these two cases. At a = 6.2 
deg, measurements were obtained at several streamwise 
stations along the airfoil and in the near wake. These data 
reveal a rapid decay in the production rate of turbulence 
kinetic energy with streamwise distance. Associated w-ith this 
rapid decay in production rate is a decay in the turbulence 
kinetic energy and Reynolds shear stress. In the near wake, 
the streamwise fluctuations decay more rapidly, w'ith a near 
isotropic condition being established at x/c- 1 . 33 . As ex- 
pected, the mixing length increases in the near wake as the 
flow establishes a wake character. The mean velocity 
measurements suggest a free-shear-layer type flow down- 
stream of separation with very rapid closure of the wake near 
the trailing edge. 

The results of this study suggest that it is at or very near the 
separation point that improved turbulence modeling is 
needed. Downstream of the separation point, a zero equation 
model would appear adequate. An accurate accounting of the 
Reynolds normal stresses at and near the separation point 
may be what is needed to predict the initial perturbations to 
the flow. 
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